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Abstract. The (J, K) = (1,1), (2,2), (3,3), and (4,4) in- 
version lines of ammonia (NH3) have been detected to- 
ward the nuclear 40" sized bar of the nearby spiral galaxy 
Maffei 2. The relative intensities of the ammonia lines are 
characterized by a rotational temperature of 85 K. This 
is higher than rotational temperatures measured toward 
IC 342 and most Galactic Center clouds, implying kinetic 
temperatures ^100 K. Since the kinetic temperature of the 
gas is larger than that of the dust, NH 3 is tracing a partic- 
ularly dense warm gas component that is heated by young 
massive stars, cloud-cloud collisions, or ion-slip heating in 
the nuclear starburst. The gas north of the nucleus (Vlsr 
= ~80kms _1 ) is more highly excited than the gas further 
south (+6kms~ 1 ). This asymmetry might be related to 
pronounced morphological distortions that are observed 
in the north-eastern part of the galaxy. 



Key words: Galaxies: individual (Maffei 2) - Galaxies: 
ISM - Galaxies: radio lines 



1. Introduction 

For the past decade, it has been possible to study the 
excitation of the molecular gas in external galaxies (e.g. 
Mauersberger & Henkel 1989; Mauersberger et al. 1990; 
Jackson et al. 1995; Huttemeister et al. 1997; Paglione 
et al. 1997; Heikkila et al. 1999; Mao et al. 2000). While 
the interpretation of CO emission often requires the ap- 
plication of complex models that incorporate temperature 
gradients at the surfaces of clouds, species like CS, HCN, 
or HC3N trace higher density gas closer to the cloud cores 
and may thus be properly modeled with the standard 
'Large Velocity Gradient' approach. Observed line inten- 
sity ratios of such high density tracers depend mostly on 
the number density n(H2), but also on the kinetic temper- 
ature Tkin of the gas. It is virtually impossible to disentan- 
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gle these parameters, making observations of a molecule 
tracing exclusively temperature highly desirable. 

Ammonia (NH3) is such a molecule. With a sin- 
gle telescope-receiver configuration at 18-26 GHz, a large 
number of NH 3 inversion transitions can be measured, cov- 
ering an enormous range of molecular excitation levels. 
These include the (J,K) — (1,1) to (4,4) inversion lines 
at 23, 65, 123, and 199 K above the ground state. Extra- 
galactic ammonia was first detected in the (J, K) = (1,1) 
line toward the nearby face-on spiral IC 342 (Martin & Ho 
1979). This was followed by detections of the (2,2), (3,3), 
and (4,4) lines (Ho et al. 1982, 1990; Martin & Ho 1986) 
that permitted a first estimate of the kinetic temperature 
of the dense nuclear gas in an external galaxy. 

In view of the four ammonia lines observed toward 
IC342 (the (6,6) line was also tentatively detected) it is 
surprising how little is known about ammonia in other 
external galaxies. Motivated by our unpublished detec- 
tions of (J,K) = (1,1) and (2,2) emission in 1994 and by 
recent improvements in receiver sensitivity, baseline sta- 
bility, and accessible bandwidth, we present a multilevel 
ammonia study of Maffei 2, an optically obscured barred 
spiral galaxy that is located behind the plane of the Milky 
Way at a distance of ^2.5 Mpc (e.g. Karachentsev et al. 
1997). High lower limits to the kinetic temperature of its 
dense nuclear gas are derived. 



2. Observations 

All data presented here were taken in March 2000, us- 
ing the Effelsberg 100-m telescope of the MPIfR equipped 
with a dual channel K-band HEMT receiver. The system 
temperature was ^220 K on a main beam brightness tem- 
perature scale; the beam size was ~40". The data were 
recorded using an autocorrelator with 8x256 channels and 
bandwidths of 80 MHz. The eight backends were config- 
ured in two groups of four, sampling data from both linear 
polarizations. Frequency shifts between the four backends 
representing a given linear polarization were adjusted in 
such a way that the (J, K) = (1,1) to (4,4) lines of ammo- 
nia could be measured simultaneously. 
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Fig. 1. NH 3 spectra towards Maffei 2 (a 20 oo = 2 h 41 m 
55fl, <S 2 ooo = +59° 36' 12") after subtraction of lin- 
ear baselines. Channel spacings are 3.9, 7.8, 3.9, and 
15.7kms -1 for the (J, K) = (1,1) to (4,4) lines at 
23694.496, 23722.631, 23870.130, and 24139.417MHz, re- 
spectively. The NH3 (4,4) spectrum, with the temperature 
scale multiplied by a factor of 2.5, is shown with a tenta- 
tive gaussian decomposition. 

The measurements were carried out in a dual beam 
switching mode (switching frequency 1 Hz) with a beam 
throw of 121" in azimuth. Calibration was obtained by 
measurements of the radio continuum of W3(OH) and 
3C 286 (for fluxes, see Baars et al. 1977; Mauersberger 
et al. 1988; Ott et al. 1994) and NH 3 spectra of Orion- 
KL and IC342 (Martin & Ho 1986; Hermsen et al. 1988). 
We estimate an absolute calibration accuracy of ±10%. 
Relative calibration should be even more accurate, at the 
order of ±5%. Pointing errors as well as weather or eleva- 
tion dependent gain variations affect all lines in a similar 
way; the calibration signal injected by a noise diode shows 
little frequency dependence and remains stable (within at 
least 5%) over the observed frequency interval. Spectral 
intensities were converted to a main beam brightness tem- 
perature (T mb ) scale. Pointing was checked every hour on 
nearby continuum sources and was found to be stable to 
within 5-10". 



Table 1. Integrated line intensities (J T m bdi>), Local 
Standard of Rest velocities (Vlsr), full width to half 
power linewidths (AVi/ 2 ), and upper state column den- 
sities (2-N u (J,K) ~ N(J,K)) for the (1,1) to (4,4) lines 
of ammonia towards Maffei 2. Line parameters were ob- 
tained from gaussian fits to the data. Errors in Cols. 3 
and 4 are standard deviations from the fit; Cols. 2 and 5 
also include a 10% uncertainty in absolute calibration (see 
Sect. 2). 



Transition 
(J,K) 


J' T mb dv 
Kkms" 1 


Vlsr 
kms 


AV 1/2 


Nn(J,K) 

10 12 cm" 2 


(1,1) 


0.527±0.080 


-78±2 


43±6 


3.44±0.52 


(1,1) 


1.024±0.121 


+6±2 


70±5 


6.69±0.78 


(2,2) 


0.694±0.122 


-78±6 


83±12 


3.40±0.60 


(2,2) 


0.782±0.131 


+5±4 


66±10 


3.83±0.63 


(3,3) 


1.127±0.131 


-82±2 


62±4 


4.88±0.57 


(3,3) 


1.110±0.128 


+8±2 


56±4 


4.80±0.55 


(4,4) 


0.882±0.215 


-57±20 


174±40 


3.55±0.87 



3. Results 

Fig. [j] shows the measured spectra toward the central 
0.5 kpc, displaying a velocity range between -530 and 
+350 kms -1 . Line parameters are given in Table |. The 
(1,1) to (3,3) lines are clearly detected; the (4,4) line is ob- 
served at a 4.5cr level. The (1,1) to (3,3) profiles show two 
distinct velocity components, separated by ~85kms -1 . In 
the (1,1) and (2,2) lines, the higher velocity component is 
stronger; in the (3,3) line, however, both components show 
similar line strength. In the (4,4) transition it is difficult 
to resolve the two features but the line center is closer 
to the velocity of the -80 kms -1 component (Table 1). 
This is also suggested by a tentative two component fit to 
the (4,4) line shown in Fig.|l], where 60-65% of the emis- 
sion seems to arise from the -80kms -1 component. We 
conclude that there is a trend with the -80kms _1 compo- 
nent rising in relative strength with excitation above the 
ground state. 

4. Determination of rotational temperatures 

With hi//k ~ 1.14K and expected excitation tempera- 
tures kJIOK across an inversion doublet we obtain, in the 
optically thin case, beam averaged column densities of 

for an individual inversion state. The column density N, 
the frequency v, and the integral are in units of cm -2 , 
GHz, and Kkms -1 , respectively. Calculated column den- 
sities are displayed in the last column of Table 1. For our 
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observed metastable (J = K) inversion transitions, the 
rotational temperature T rot between the states (J, J) and 
(J 7 , J') can then be obtained from 

NjJ, J) = g op (J') 2J' + 1 f- AE \ 
N(J',J') g op (J) 2J+1 eXp \kT mt )' 

where g ap = 1 for para-ammonia, i.e. for (J,K) = (1,1), 
(2,2), and (4,4), and g op — 2 for ortho-ammonia, i.e. for 
the (3,3) line. AE is the energy difference between the 
states involved. 

Rotation diagrams of the ammonia emission in the 
(1,1) to (4,4) lines are shown in Fig.|. The (3,3) line, al- 
though belonging to a different ammonia species, fits well 
in the general trend of decreasing normalized column den- 
sity with energy above the ground state. The slope of the 
lines fitting the normalized column densities is -loge/T rot . 
A weighted fit to the total ammonia emission (solid line; 
relative weights reflect signal-to-noise ratios) yields a rota- 
tional temperature of 85±15K (la error). Excluding the 
(4,4) line from the fit gives 65±10K. For the two veloc- 
ity components, only including the (1,1) to (3,3) lines, 
T rot = 85±5 (Vlsr = -SOkms" 1 ) and 55±10K (F LSR = 
+6kms _1 ). Inclusion of the highly tentative decomposi- 
tion of the (J, K) = (4,4) line (Fig.0) yields T rot ~ 115 
and 65 K. As sometimes observed in Galactic sources, the 
inclusion of higher excited inversion levels leads to higher 
rotational temperature estimates. This is indicated by the 
decreasing (negative) slope of the line connecting the mea- 
sured points. 

5. Discussion 

5.1. Interpretation of rotational temperatures 

The J = K ammonia inversion doublets are called 
metastable because they decay radiatively via slow AK 
= ±3 transitions. As a consequence, the relative popula- 
tion of these states is strongly affected by collisions and 
thus reflects (e.g. Walmsley & Ungerechts 1983) the ki- 
netic temperature of the gas (T rot < Tkin)- 

One way to explain rising rotational temperatures with 
increasing energy above the ground state is the pres- 
ence of kinetic temperature gradients. Alternatively, op- 
tically thick lines or subthermal excitation can produce 
the same effect (Goldsmith & Langer 1999). Accounting 
for these latter two effects, radiative transfer calculations 
(e.g. Walmsley & Ungerechts 1983; Schilke 1989) predict 
that rotational temperatures approach the kinetic temper- 
ature with rising energy of the analysed metastable inver- 
sion doublets. This yields rising T rot with rising E/k. It 
is therefore likely, based on the corrections of the models, 
that the true kinetic temperature is >100K. Tki n should 
be larger for the -80 than for the +6kms~ 1 velocity com- 
ponent. 
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Fig. 2. A rotation diagram of metastable ammonia tran- 
sitions towards Maffei 2, assuming optically thin emission 
and excitation temperatures 3>2.7K for each inversion 
doublet. The column densities N from Table 1 refer to 
the upper state of the respective inversion doublet. Filled 
squares: Total emission (dashed line connects the points); 
empty squares: -80kms _1 component; empty triangles: 
+6kms _1 component. Straight lines provide best fits to 
the total emission in the (J, K) — (1,1) - (4,4) lines (solid) 
and to the (1,1) - (3,3) emission lines for the -80 and 
+6kms _1 components (dashed and dotted, respectively). 
The decomposition of the (4,4) line into two velocity com- 
ponents, based on the fit shown in Fig. [l], is at best only 
marginally significant. Errors for the individual (1,1) to 
(3,3) components are ~40% larger than those for the to- 
tal emission. 



5.2. A comparison with IC 342 and Galactic Center 
clouds 

NH3 integrated line temperatures are larger towards Maf- 
fei 2 than towards IC 342, also located at £>~2.5 Mpc 
(Karachentsev et al. 1997). Emission from the inner 
500 pc of the Milky Way would show similar intensities 
if observed from this distance (e.g. Hiittemeister et al. 
1993a,b). 

From the (J, K) = (1,1) to (4,4) lines, T rot ~85K in 
Maffei 2, but only 50±10K in IC 342 (Martin & Ho 1986). 
Toward peak positions of clouds near the center of the 
Galaxy, (1,1) lines tend to be optically thick, while more 
highly excited inversion lines show no significant optical 
depth effects (Hiittemeister et al. 1993b). There are two 
main temperature components, one at Tki n ~ 25 and one 
at ~200K. These, if present in Maffei 2 as well, could 
be responsible for the increasing rotational temperatures 
with E/k (Fig. ^J). Such gradients are also seen in IC342, 
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if the (6,6) line is as strong as suggested by the spectrum 
presented by Martin & Ho (1986). 

5. 3. The nuclear region of Maffei 2 

What kind of molecular gas is seen in ammonia? Critical 
densities, where collisional de-excitation matches sponta- 
neous emission rates, are n(fi2) ~ 3T0 3 cm -3 (e.g. Green 
1980; Ho & Townes 1983). Our beam size (40") matches 
the central bar with its size of 43" x 9" in 12 CO J = 1-0 
(Ishiguro et al. 1989) and 40"xl2" in 13 CO J = 1-0 (Hurt 
& Turner 1991). The presence of bars provides a viable 
mechanism to transfer interstellar gas into the nuclear re- 
gion of a galaxy, leading to the formation of dense molecu- 
lar clouds and triggering bursts of star formation. As long 
as the gas is not yet fully ionized or swept away by stel- 
lar winds, it will coexist with the newly formed stars and 
will be heated by their radiation, by an enhanced flux of 
cosmic rays, by cloud-cloud collisions, and by the dissipa- 
tion of turbulent energy. Both Maffei 2 and IC 342 show 
prominent nuclear bars and strong CO, infrared, and ra- 
dio continuum emission (e.g. Lo et al. 1984; Ishiguro et al. 
1989; Ishizuki et al. 1990; Hurt & Turner 1991; Turner & 
Ho 1994; Meier et al. 2000). 

Assuming a grain emissivity law of A -1 , dust tem- 
peratures deduced from the Kuiper Airborne Observa- 
tory and, for a larger region, from IRAS data (Rickard 
& Harvey 1983; IRAS 1989) are T d = 33-47 K. Ignor- 
ing non-metastable inversion doublets, Table 1 implies 
that iV(NH 3 ) is on the order of 10 14 cm -2 . For H 2 , with 
£>=2.5Mpc, AT(H 2 )~10 22 cm- 2 (Hurt & Turner 1991). 
This yields a relative abundance of A"(NH 3 )^10~ 8 . With 
T ro t(NH3) ^> Td and an unusually small NH3 abundance 
(e.g. Benson & Myers 1983; Brown et al. 1988), ammo- 
nia should trace a gas component that is different from 
that seen in CO (for an analysis of CO emission in a nu- 
clear starburst environment, see Mao et al. 2000). Cloud- 
cloud collisions, young massive stars, or ion-slip heating 
(Hiittemeister et al. 1993a, b) may cause the high rota- 
tional temperatures. SiO, an exclusive tracer of high tem- 
perature gas, was detected in Maffei 2 (Sage & Ziurys 
1995) and is likely arising from the same gas component. 
Detailed studies of NH3 and SiO in the Galactic Center re- 
gion (e.g. Hiittemeister et al. 1993a, b, 1998) indicate that 
the warm gas not directly associated with massive stars is 
only moderately dense (n(H2)~5T0 3 cm~ 3 ). 

From a comparison with the interferometric maps 
of Ishiguro ct al. (1989) and Hurt & Turner (1991), 
the ammonia component at +6 kms must arise SW, the 
-SOkms^ 1 component NE of the center of the galaxy. Why 
is the rotation diagram (Fig.|2j) indicating higher excita- 
tion NE of the nucleus? The bar in Maffei 2 may be caused 
by a merger with a satellite galaxy in the NE that severely 
disrupts this side of the galaxy (Hurt et al. 1996). This 
might have led to differences in inflow rates, star forming 
rates, and gas heating at the north-eastern and south- 



western sides of the nucleus. While our observed small 
scale (~500pc) nuclear asymmetry in rotational temper- 
ature may be consistent with the large scale asymmetry 
observed in Hi, the proposed connection between these 
phenomena remains speculative. More data on the puta- 
tive interacting dwarf galaxy, detailed simulations of the 
interaction, and models of the evolution of the nuclear bar 
are needed to verify or to reject the proposed scenario and 
to explain the recent past of Maffei 2. 

Acknowledgements. We wish to thank J.L. Turner for critically 
reading the manuscript. 
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